Decay of scalar turbulence revisited 
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We demonstrate that at long times the rate of passive scalar decay in a turbulent, or simply 
chaotic, flow is dominated by regions (in real space or in inverse space) where mixing is less efficient . 
We examine two situations. The first is of a spatially homogeneous stationary turbulent flow with 
both viscous and inertial scales present, ft is shown that at large times scalar fluctuations decay 
algebraically in time at all spatial scales (particularly in the viscous range, where the velocity is 
smooth). The second example explains chaotic stationary flow in a disk/pipe. The boundary region 
of the flow controls the long-time decay, which is algebraic at some transient times, but becomes 
exponential, with the decay rate dependent on the scalar diffusion coefficient, at longer times. 

PACS numbers: 47.27.Qb 



Study of advection of a steadily supplied scalar in 
smooth flow, pioneered by Batchelor [1] for the case of a 
slowly changing strain and extended by Kraichnan to the 
opposite limit of short-correlated flow [2] , has developed 
into a universal theory, applicable to any kind of statis- 
tics and temporal correlations of the smooth chaotic flow 
[3, 4]. (See also the reviews [5, 6].) Extension of this 
theory to the decay problem [7, 8] suggested a fast (ex- 
ponential) decay of the passive scalar in the case of spa- 
tially smooth (that is approximated by linear velocity 
profiles at relevant scales) and unbounded random flow. 
(We will call this case the "pure Batchelor" one.) In par- 
allel, laboratory experiments [9, 10, 11, 12] and numeri- 
cal simulations [13, 14, 15, 16] were conducted to verify 
the Batchelor-Kraichnan theory. However, a compari- 
son between theory, experiment and numerical simula- 
tions was not conclusive. Different experiments and sim- 
ulations showed different results, inconsistent with each 
other and with the theory. In this letter we explain a pos- 
sible source of the discrepancy. Our major point is that 
the Batchelor-Kraichnan theory does not apply to the 
late-time stage of the scalar decay measured in the de- 
cay laboratory and numerical experiments [9, 10, 12, 13]. 
This is because the long-time decay is dominated by the 
regions of the flow, where the mixing is not as efficient 
as in the pure Batchelor case. 

Before constructing a quantitative theory, let us draw 
a qualitative physical picture. We start from the first ex- 
ample of a multi-scale turbulent flow with the following 
hierarchy of scales assumed: Ly ^ rj ^ r^, where Ly 
is the energy containing scale of the flow, rj is the Kol- 
mogorov (viscous) scale and is the diffusive scale of the 
scalar. It is known that the passive scalar mixing in the 
inertial range of turbulent flow is slower at larger scales, 
since the turnover time grows with scale. This explains 
the algebraic-in-time decay of the scalar for the inertial 
range of scales, Ly ^ r ^ -q [17, 18]. If we now take 
into account that turbulent velocity is smooth at scales 
smaller than 77, it becomes important to understand what 



happens with the scalar fluctuations at these small scales. 
We show in the letter that scalar decay is accompanied 
by both upscale and downscale transport, however, the 
situation is asymmetric. The decay of scalar fluctuations 
in time is due to inertial interval physics, while the spatial 
correlations of the scalar at the smallest scales are mainly 
controlled by the viscous part of the flow. (Here, the spa- 
tial correlations of the scalar are logarithmic, like in the 
stationary case [4].) An analogous situation is realized 
if the scalar is advected by a chaotic (spatially smooth) 
flow in a box or a pipe. Advection of the scalar slows 
down in the boundary region, while mixing continues to 
be efficient in the bulk, i.e. far from the boundaries. The 
result is that the boundary region supplies the scalar into 
the bulk. Plumes (blobs) of the scalar, injected from the 
boundary domain into the box center, cascade in a fast 
way down to where diffusion smears it out. Therefore, 
the temporal decay of the scalar is mainly controlled by 
the rate of scalar injection from the boundary layer. This 
causes an essential slow down of the decay in compari- 
son with the pure Batchelor case, since in this case there 
are no stagnation (boundary) regions. The decay of the 
scalar is algebraic during the transient stage, when the 
scalar is supplied to the bulk from the boundary region 
which is wider than the diffusive boundary layer. The 
boundary region width decreases with time. Whenever 
the boundary layer width becomes of the order of the 
diffusive layer width the decrease stops and the algebraic 
temporal decay turns into an exponential one. The expo- 
nential decay rate appears to be parametrically smaller 
than in the pure Batchelor case. 

In all the cases discussed, we consider decay of a passive 
scalar, 6, described by the equation 



(1) 



where v is the flow velocity and k, is the diffusion coef- 
flcient. Our goal is, starting from the equation (1) and 
assuming some statistical properties of the velocity held 
V, to describe statistics of the passive scalar decay. In 
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order to understand the qualitative features of the pas- 
sive scalar decay, we choose (in the spirit of Kraichnan 
[19]) the simplest possible case of a short-correlated in 
time velocity, possessing Gaussian statistics. This model 
enables us to obtain analytically a time evolution of the 
passive scalar correlation functions. Two situations are 
discussed. The first is of a statistically homogeneous 
flow (infinite volume) with a multi-scale velocity which is 
smooth at small scales, and possesses a scaling behavior 
at larger scales. The Reynolds number Re = (L^/t])''/^ 
and the Prandtl number Pr = {rj/rd)^ are both assumed 
to be large. The second situation is 2d chaotic flow con- 
fined within a disk. We also discuss an extension to the 
case of a more realistic pipe fiow. The fiow is assumed 
to be smooth and obeying the standard viscous behav- 
ior in the neighborhood of the boundary: the longitu- 
dinal and transverse components of the velocities tend 
to zero quadratically and linearly (respectively) with the 
distance to the boundary. 

We begin with an analysis of the homogeneous un- 
bounded case. The statistics of a short-correlated in time 
Gaussian velocity is completely characterized by its pair 
correlation function 

{va{ti,r)vp{t2,0)) = 6{ti - t2)[VoSo,f^ - ICc^pir)] , (2) 

K^0^Lw'{r)(s^0-'-^) +^W{r)S^^, (3) 

where d is dimensionality of space. Vq is a constant 
characterizing fiuctuations of the velocity at the integral 
scale, Ly, and W{r) is a function, which determines the 
scale-dependence of velocity fluctuations. To describe 
both viscous and inertial ranges, we introduce 

W = £>rV[l + (r/r?)^] , (4) 

where D stands for the amplitude of the velocity fluc- 
tuations and 7 characterizes the velocity scaling in the 
inertial interval. The Kolmogorov scale r] separates the 
inertial and viscous intervals. The diffusive length is now 
expressed as = \J KjD. Measuring time, i, in the 
units of 1/f, one replaces D by unity in all the forth- 
coming formulas. The object we examine here is the si- 
multaneous pair correlation function of the passive scalar, 
F{t,r) = {6{t,r)0{t,0)), which is the simplest non-zero 
correlation function in the homogeneous case. Under the 
assumption of isotropy one derives from Eqs. (1-3) 

- r^-'^dr {r-^-i [W{r) + k] drFx} = XFx, (5) 

where F\(r) is the Laplace transform of F{t,r) with re- 
spect to t. The operator on the left-hand side of Eq. (5) 
is self-adjoint and non-negative (with respect to the mea- 
sure J drr'^~^ and imder condition that F{r) is smooth 
at the origin and drF{r = 0) =0). 

We assume that Jdr6 = (the Corrsin invariant 
is zero), and also, that initially the scalar is corre- 
lated at the scale ro from the viscous-convective range. 



rd ^ ro rj, i.e. F(0, r) is ss F(0, 0) at r < ro and 
decays fast enough at r > tq. Then some part of the 
initial evolution of F{t, r) is not distinguishable from the 
decay in the pure Batchelor case, described by the limit 
ry ^ oo in Eq.(4). Let us briefly describe this initial 
stage of the dacay. The solutions of Eq. (5) are power- 
like. Fa oc 7.±-\/<iV4-A-d/2^ outside the diffusive range, 
r ^ rd,. The functions are normalizable only if A > d'^/A. 
Therefore there is a gap in the spectrum (0 < A < 
are forbidden). The existence of this gap guarantees an 
exponential decay of F{t, r) with time. This feature is, 
of course, in agreement with the previous analysis of the 
pure Batchelor case [7, 8]. A complete set of functions 
at r rd is Fx = r~'^^^ cos[ak + kln{r/rd)], where k 
is a positive number related to A via A = d^/4: + k"^, 
and ttk are phases, determined by matching at r ~ r^. 
Taking into account the orthogonality condition for Fx, 
J drr'^~^Fx{r)Fx'{r) = nS{k — k')/2, one can express 
F{t,r) via F{0,r). Analysis of this expression leads to 
the following picture. The correlation function, F(t,r), 
does not change at r < r_ = ro exp(— td). At the scales 
larger than r_, F decays exponentially, oc exp[— frf"/4]. 
Note, that r_, also, marks the position of the maximum 
of the spatial derivative of F{t,r), drF{t,r). Therefore, 
r-{t) describes the front running downscale (from ro). 
The front reaches r^ at td = ln{ro/rd)/d, so that at 
t > td, F{t, r) decays exponentially at all the small scales. 
Complimentary to the front running downscale, there ex- 
ists another front running upscale. Indeed, the integral 
quantity, Ar r'^~^ F{t,r), which determines the over- 
all amount of the scalar at all the scales smaller than s, 
achieves its maximum around r_|_ = roCxp(t(i). 

At i > t,, ~ \u[r]/ro]/D the pure Batchelor description 
ceases to be valid and one does need to account for com- 
plete multiscale form of eddy-diffusivity function W{r) 
given by Eq. (4). Our further analysis is devoted to this 
general case. The spatial decay of the functions Fx in 
the multi-scale case is not as steep as in its pure Batch- 
elor one. It is convenient to change to a new field, 
Fx = r'^-'^dr[r'^/'^^x{z)], where z = 2y/Xr^/'^ /-f. The 
general solution of the resulting equation for ^x{z) is a 
linear combination of the Bessel functions, $a = J±iyiz), 
where f = 2y/d'^/4 - A/7. If A < the positive root 
for <I>A has to be chosen to satisfy the matching condi- 
tions at rd- The eigen function is normalizable for any 
positive A, so that even the smallest positive A are not 
forbidden. The lack of the gap in the spectrum of the 
multi-scale model means an algebraic in time decay of 
F{t,r), i.e. much slower decay than the one found for 
the Batchclor-Kraichnan model. We present here the ex- 
pression for the long-time, t ^ tn, asymptotic behavior 
of F(t,r): 

F(f rw/ [td + Hr/^) + {rhVh]~^~'"\ r « r+: , . 

i'(t,r)(X<^ ^_5/4-d/27j.-<i/2-H7/4^ r > r+ ^ ' 
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Here r+ = (7^)^^^^ stands for position of the front run- 
ning upscale, r+ lies in the scaling region r+ ^ rj, where 
the multi-scale model turns into the so-called Kraichnan 
model [19], W r"^'^ . (Therefore, it is not surpris- 
ing, that our result (6) is consistent at r > r+ with the 
expression for the pair correlation function derived be- 
fore for the Kraichnan model [17, 18] in the regime of 
zero Corrsin invariant.) Note, that when t — tr/ is yet 
moderate in value, the small-scale part of the asymptotic 
expression (6) is correct only at r ':$> R-{t) = rie~*'^, 
where R- is the position of the front initiated at f ~ 
and running downscale from rj. At t tr, + \n{ri/rd)/d, 
the downscale front reaches the dissipativc scale, and 
afterwards the expression (6) is correct for any scales 
r » r,;. Eq. (6) shows that just like in the pumped 
case [1, 2, 4] the scalar structure fimction is logarithmic, 
{[0{t,r)-0{t,O)]^) ~ i-2-<i/')'ln[r/rd], at < r < r? and 
large times, where, therefore, the time-dependent factor 
at the logarithm can be interpreted as a scalar flux from 
the inertial range, r > 77, downscale to the viscous range. 

To conclude, the existence of the inertial region in ve- 
locity affects drastically the spatio-temporal distribution 
of the scalar in the viscous range (where the velocity can 
be approximated by linear profiles). The inertial inter- 
val, where the scalar is mainly concentrated, serves as a 
kind of reservoir for the smaller scales. This explains why 
the decay of scalar correlation function in the viscous is 
much slower (algebraic) than in the pure Batchelor case 
(when it would be exponential). 

Now wc proceed to the situation of spatially bounded 
flows. Aiming to demonstrate the main qualitative fea- 
tures of the confined geometry, that is strong sensitivity 
of the scalar decay to the peripheral (close to the bound- 
ary) part of the flow, we examine a model case of 2d 
chaotic (i.e. consisting of only few spatial harmonics) 
flow inside a disk. (It is, actually, clear that the picture 
of temporal decay and spatial distribution of the scalar, 
described below, is imiversal, i.e. it applies to a typi- 
cal chaotic; flow confined in a close box, particularly in 
'id.) Incompressible flow in 2d can be characterized by 
a stream function ip, then the radial and azimutal com- 
ponents of the velocity are Vr = —d^tp/r and v^p = drip- 
Our model of ip is 

i> = ~r'^U{r) sin(2(^) + ^r^U{r) cos(2(/?) , (7) 

{^i{ti)^j{t2)) =2D6ijSiti-t2), (8) 

where U (r) is a function of r, finite at the origin, 
U{0) = 1, and becoming zero, together with its first 
derivative, at the disk boundary, r = 1, and ^1,2 (t) 
arc zero mean short correlated random Gaussian func- 
tions. The value of the passive scalar, 9, averaged over 
the statistics of ^, {6{t,r)), is the object of our inter- 
est here. One examines how the average concentration 
of the scalar evolves with time at different locations r 



within the disk. The short-correlated feature of the ve- 
locity field allows a closed description for (6*) in terms of 
a partial differential equation. Considering the spheri- 
cally symmetric part of (9) only (asymptotically, at large 
times only this uarp/iz-independent part remains essen- 
tial), and passing from (9) and r to T and q, respec- 
tively, where (9) = r~^dr[r'^T{t, q)] and q = — Inr, one 
finds that the Laplace transform of T with respect to t 
(measured, again, in units) satisfies 

{U^ + Ke-29)(a2TA - 2a,TA) + ATa = 0. (9) 

The diffusion term is important only at the center of the 
disk and near the boundaries. Analysis of Eq. (9) is 
straightforward but bulky. Below we will present only 
selected details of the analysis, aiming to describe the 
general picture of the phenomenon. (The complete ac- 
count for the derivation details will be published else- 
where [20].) 

In the bounded flow the decay of passive scalar splits 
into three distinct stages. The major effect dominating 
the first stage (just as in the pure Batchelor case, ex- 
plained by Eq. (9) with {/ — > 1) is formation of elon- 
gated structures (stripes) of the scalar in the bulk region 
of the flow. The stripes are getting thinner with time, 
i.e. inhomogeneities of smaller and smaller scales are 
produced. Once the width of the stripes decreases down 
to the dissipative scale, ~ ^/k, the stripes are smeared 
out by diffusion. The stretching-contraction process is 
exponential in time, so that the initial stage when the 
stripes are formed lasts for n ~ ln[l/rd]. By the end 
of this first stage the scalar is exhausted in the cen- 
tral region of the flow. The stretching rate, however, 
is smaller in the peripheral domain than in the bulk. 
Thus the scalar literally remains longer in the periph- 
eral domain. This deflnes the second, transient, stage. 
At K^/"* < g « 1 - r < X^/* < 1, the solution of Eq. 
(9) is Tx oc qsm^ax + \/Xq/U^. In the other asym- 

potitc region atq^ A^/* one can drop the second deriva- 
tive (with respect to q) and diffusive terms in Eq. (9). 
One flnds that, 9^ In ~ 1 at g ~ A^^^. Matching 
those two asymptotic regions at g A^/* we come to the 
ax = condition. Once the asymptotic form of the eigcn- 
function is known, it is straightforward to restore the 
behavior of T{t,q). One finds that T{t,q) (and, there- 
fore, (9)) is concentrated in the S{t) ^ l/\/t-small (and 
shrinking with time) vicinity of the boundary. In the do- 
main outside of the shrinking layer the decay is algebraic, 
(9) (X t^'^/^g^^. The second stage lasts for T2 ~ k^-'^/^, 
i.e. until S{t) shrinks down to ru = k^^^, which is the 
width of the diffusive boundary layer. Account for diffu- 
sivity in the q ^ 1 analysis gives yet another matching 
condition for T^, at r ~ r^, resulting in formation of a 
discrete spectrum for A. The smallest A (and therefore 
the value of the level spacing in the discrete spectrum) 
is estimated by y/R,. (Using the periodic character of 
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the sin-function, one finds that the n-th level eigenvalue 
A„ is estimated by ~ n^^/H at large n.) Therefore, at 
i ^ T2, {&) decays exponentially, oc exp{—t/td), where 
td ~ xhe smallest A eigen-hmction is localized 

at g ~ ru- Thus the third (final) stage of evolution is 
characterized by the majority of the scalar remaining in 
the diffusive boundary layer. 

To conclude, the temporal behavior of the passive 
scalar correlations, predicted by the bounded flow the- 
ory, is complicated. It involves two different exponential 
regimes separated by an algebraic one, so that a spe- 
cial accuracy is required in order to quantify the theory 
against various experiments in the field. 

The bounded flow theory can be applied to the exper- 
iment of Groisman and Steinberg [12], where the pas- 
sive scalar is advected through a pipe by dilute-polymcr- 
solution flow. The polymer-related elastic instability 
makes the flow chaotic. The mean flow along the pipe 
is also essential, so that according to the standard Tay- 
lor hypothesis, measurements of scalar concentration at 
different positions along the pipe are interpreted as cor- 
respondent to different times in an artificial decay prob- 
lem. Exponential decay of the scalar was reported in [12] 
at long times (long pipes). The theoretical picture ex- 
plaining the experiment is similar to the one proposed 
above for the disk. The width of the diffusive boundary 
layer near the pipe boundaries is estimated by R-Pe~^^'^, 
where R is the pipe radius and Pe is the Pcclet number, 
Pe = R^a/K and cr is a typical value of the velocity gra- 
dient (Pe is ~ 10^ in the conditions of [12]). The charac- 
teristic decay time (fixed by the diffusive boundary layer 
width) is estimated by td ~ a^^Pe^^^. However, one 
should be careful when converting this time to the decre- 
ment of the passive scalar decay along the direction of 
the fiow (pipe). The average velocity of the fiow tends 
to zero near the boundary of the pipe, making the ad- 
vection near the boundary less efficient than in the bulk. 
Assuming a linear profile of the average velocity near the 
boundary, one gets the following estimate for the law of 
scalar fluctuation decay along the pipe, ^ cxp(— 72/77,0), 
where z is the coordinate in the direction of the pipe, 
uo is the average velocity at the center of the pipe, and 
7 ^ (T Pe~^^^. The factor Pe^^^^ here is a manifestation 
of a slow down of the passive scalar decay in comparison 
with the pure Batchelor case. The same factor charac- 
terizes decay of higher order correlation functions of the 
scalar as well. For example, one can consider the passive 
scalar pair correlation function F{r) near the center of 
the disk. One finds that F{r) cx r~". where a is small 
and can be estimated as a ~ Pe"^/''. (Note, that this 
behavior is almost indistinguishable from the logarith- 
mic one, found for the stationary case [1, 2, 4].) This 
estimations agree with the experimental data of [12]. 

Let us now briefly discuss two other available numerical 
and experimental results. Pierrehumbert [13] reported 



exponential decay of the scalar correlations in his numer- 
ical experiment with chaotic map velocity in a periodic 
box. This observation is consistent with our results, since 
the periodic boimdary conditions arc less restrictive than 
the zero condition for velocity at the box boundary, and, 
therefore, should leave a finite gap in the A-spectrum even 
in the limit k 0. In the experiment of JuUien, Cas- 
tiglione and Tabeling on 2d stationary flow steered by 
magnets in a flnite volume beneath the 2d layer [11], the 
decay rate of the scalar seems slower than exponential 
(Fig. 2 of [11]). This observation is in agreement with 
the absence of the gap in the A-spectrum (at k — > 0) we 
found for the finite box model. 

The brevity of this letter does not allow us to discuss ef- 
fects of intermittency (which manifests in high-order mo- 
ments of the scalar), of higher order angular harmonics 
(in the disk or pipe geometry), and details of the spatial 
distribution of the scalar in a variety of other inhomoge- 
neous cases (e.g. periodic flow). This detailed discussion 
is postponed for a longer paper to be published elsewhere 
[20] (where we will also present results of direct numerical 
simulations of the passive scalar decay for these cases) . 
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cussions. 



[1] G. K. Batchelor, JFM 5, 113 (1959). 

[2] R. H. Kraichnan, Phys. Fluids. 10, 1417 (1967); JFM 47, 

525 (1971); Ibid 67, 155 (1975). 
[3] B.I. Shraiman, E.D. Siggia, Phys. Rev. E 49, 2912 (1994). 
[4] M. Chertkov, G. Falkovich, I. Kolokolov, and V. Lebedev, 

Phys. Rev. E 51, 5609 (1995). 
[5] B. I. Shraiman and E. D. Siggia, Nature 405, 639 (2000). 
[6] G. Falkovich, K. Gaw§dzki, and M. Vergassola, Rev. 

Mod. Phys. 73, 913 (2001). 
[7] D. T. Son, Phys. Rev. E 59, R3811 (1999). 
[8] E. Balkovsky, A. Fouxon, Phys. Rev. E 60, 4164 (1999). 
[9] P. L. Miller and P. E. Diniotakis, JFM 308, 129 (1996). 
[10] B. S. WiUiams, D. Marteau, and J. P. GoUub, Phys. Flu- 
ids 9, 2061 (1997). 
[11] M.-C. JulUen, P. Castiglione, and P. TabeUng, Phys. Rev. 

Lett. 85, 3636 (2000). 
[12] A. Groisman and V. Steinberg, Nature 410, 905 (2001). 
[13] R. T. Pierrehumbert, Chaos. Sol. Fract. 4, 1091 (1993). 
[14] M. Holzer and E. Siggia, Phys. Fhiids 6, 1820 (1994). 
[15] T. M. Antonsen, F. F. Fan, and E. Ott, Phys. Rev. Lett 

75, 1751 (1995). 
[16] D. Bogucki, J. A. Domaradzki, and P. K. Yeung, JFM 

343, 111 (1997). 
[17] G. Eyink and J. Xin, J. Stat. Phys. 100, 679 (2000). 
[18] M. Chaves, G. Eyink, U. Frisch, and M. Vergassola, Phys. 

Rev. Lett. 86, 2305 (2001). 
[19] R. H. Kraichnan, Phys. Rev. Lett. 72, 1016 (1994). 
[20] M. Chertkov, V. Lebedev, and Z. Toroczkai, in prepara- 
tion. 



